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1  INTRODUCTION. 

1.1  THE  SPIN  PROBLEM. 


The  stall  departure  of  general  aviation  aircraft  is  a  major 
contributor  to  aircraft  accidents.  The  desire  for  designing  new 
aircraft  to  be  spin  resistant  has  been  revived  in  recent  years, 
in  order  to  reduce  the  chance  of  inadvertent  spin  entry  is 
greatly  reduced.  Another  goal  is  to  provide  pilot  warning  when 
the  conditions  are  favorable  for  spin  entry,  and  provide  it  early 
enough  to  enable  him  to  take  corrective  action  before  aircraft- 
control  or  excessive  altitude  are  lost. 

The  work  described  herein  deals  primarily  with  the  first  goal 
with  two  different  but  related  efforts. 

1.  Development  of  a  real-time  visual  simulator  with 
interactive  controls  to  enable  evaluation  of  the 
computer  code  used,  as  well  as  to  evaluate  aircraft 
configurations  being  designed. 

2 .  Development  of  an  aerodynamic  prediction  method  for  wing 
planforms  to  enable  the  evaluation  of  the  stall/spin 
behavior  computationally. 

1.2  CRITERIA  FOR  THE  DEVELOPMENT  WORK . 

The  techniques  developed  are  simple  enough  that  advanced  concept 
teams  with  limited  computing  facilities  and  budgets  can  take 
advantage  of  them.  This  means  that  the  visual  simulation  should 
be  performed  on  a  standard  high  performance  graphics  workstation, 
either  stand-alone  or  in  conjunction  with  a  host  computer.  Near 
real-time  simulation  was  also  necessary  so  that  pilot  inputs  are 
easy  to  apply. 

For  the  aerodynamic  prediction  code,  the  method  should  be  based 
on  a  well-known  technique  of  solving  air  flow  problems  and 
requiring  modest  computing  resources.  Again,  this  is  keeping  the 
configuration  design  teams  in  mind. 

2_j_  THE  SIMULATION. 

2 . 1  EQUATIONS  OF  MOTION. 

In  the  stall/spin  entry  regime,  the  full  non-linear  coupled 
equations  of  motion  must  be  solved  as  a  marching  problem.  The 
aerodynamic  forces  and  moments  are  functions  of  the  motion 
variables,  namely  angles  of  attack  and  sidesliD,  velocity  and 
angular  velocity.  For  this  reason,  the  aerodynamic  lift  and  drag 
are  integrated  at  each  timestep  over  the  wing  span. 

Details  of  these  computations  are  in  Appendix  A,  a  reproduction 
of  a  raster's  Thesis  by  J.  Ritter  on  the  subject. 
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2.2  VISUAL  SIMULATION. 

The  results  of  the  numerical  solution  of  the  equations  of  motion 
at  each  time  step  are  the  three  spatial  coordinates  and  the  three 
orientation  angles  of  the  aircraft. 

These  six  numbers  are  used  to  perform  the  necessary  geometric 
transf ormation  for  the  scene  to  be  displayed  on  the  workstation 
screen.  Three  options  are  available  and  discussed  below. 

2.2.1  View  From  the  Cockpit . 

A  scene  of  the  surrounding  environment  was  developed  to  give  the 
operator  a  feeling  for  the  rotations  and  altitude  at  the  current 
time.  The  environment  consists  of  a  flat  terrain  with  a  river 
and  valley.  Also  included  is  a  sea,  two  ranges  of  mountains,  a 
sky  and  some  scattered  clouds.  This  is  a  simple  enough  scene  to 
expedite  the  processing,  yet  give  a  good  feel  for  the  aircraft 
motion . 

An  instrument  panel  comprised  of  speed,  altitude  and  rate  of 
climb  indicators,  an  artificial  horizon  and  a  compass  were  added. 
Also  included  was  a  grid  to  indicate  the  position  of  the  three 
controls:  rudder,  elevator  and  aileron. 

It  was  learned  that  including  the  instrument  panel  display 
required  an  excessive  amount  of  processing  time,  thus  making  it 
impossible  to  achieve  real-time  simulation.  It  was  therefore 
removed  temporarily,  until  code  optimization  can  be  performed. 

2.2.2  View  From  The  Ground. 

A  view  of  the  environment  scene  plus  the  aircraft  from  a  fixed 
point  on  the  ground  (zero  altitude)  can  be  achieved  in  this 
option.  A  hidden  surface  algorithm  was  developed  to  display  the 
aircraft  realistically.  The  method  is  to  generate  an  aircraft- 
fixed  coordinate  system,  and  at  each  time  step,  determine  which 
octant  of  the  coordinate  system  the  view  point  occurs.  One  of 
eight  possible  sequences  for  drawing  the  panels  of  the  aircraft 
is  then  performed,  depending  on  which  octant  is  selected.  A  zoom 
feature  is  also  employed. 

2.2.3  View  From  A  Flying  Observer. 

Because  of  limits  on  the  zoom  angle  in  the  above  option,  the 
aircraft  appears  very  small  when  at  high  altitudes.  A  moving 
viewpoint  was  provided  as  a  third  option  to  simulate  the  view 
from  an  observation  aircraft  flying  parallel  to  the  object 
aircraft,  except  it  does  not  spin  or  rotate.  This  viewpoint  is 
displaced  from  the  object  aircraft  by  a  fixed  vectorial  length. 
In  addition  to  being  closer,  thus  viewing  the  aircraft  better, 
the  effect  of  change  in  the  line  of  sight,  which  could  be 
misleading,  is  eliminated.  This  also  simulates  the  observers 
view  in  a  spin  tunnel  test.  All  the  other  features  of  Section 
2.2.2  are  included  here. 
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2 . 3  HARDWARE . 

Even  though  it  was  originally  intended  to  perform  the  solution  of 
the  equations  of  motion  on  a  host  computer,  limitations  on  the 
data  transfer  rates  of  the  Ethernet  made  it  difficult  to  modify 
the  program  structure  to  take  advantage  of  certain  features  of 
the  workstation.  Alternatives  include  Direct  Memory  Access  (DMA) 
devices,  which  limit  the  physical  distance  between  the  host  and 
the  terminal. 

The  ideal  solution  may  be  to  use  a  super  workstation  with  higher 
computation  speed  and  to  obtain  advanced  training  for  taking  full 
advantage  of  the  system  capabilities.  A  full-scale  Learjet 
cockpit  simulator  was  acquired  from  Learjet  Company  recently. 
This  will  be  fitted  with  visual  display  and  will  employ  the 
cockpit  controls  in  future  work. 

3 .  COMPUTATIONAL  EVALUATION  OF  WING  PLANFORMS . 

Most  computer  codes  available  today  are  either  limited  to  low 
angle  of  attack,  small  perturbation  aerodynamic  analysis,  or 
require  considerable  modeling  and  computing  power.  An  example  of 
the  latter  is  the  Navier-Stokes  codes  which  run  on  super 
computers  and  are  difficult  to  use.  These  codes  have  not  yet 
reached  the  level  of  reliability  which  would  permit  preliminary 
designers  to  use  them. 

This  part  of  the  current  effort  provides  a  method  which  uses  the 
familiar,  efficient,  and  reliable  panel  codes  with  minor 
modifications.  This  technique  and  its  application  to  a  wing  are 
detailed  in  Appendix  B. 

4 .  FUTURE  WORK. 

Building  on  the  experience  of  current  work,  the  following 
objectives  are  identified  to  provide  industry  and  government  with 
a  research  and  development  facility  which  deals  with  pre¬ 
prototype  and  pre-wind  tunnel  study  of  the  spin  entry  problems. 
This  will  enable  configuration  designers  to  evaluate  their 
concepts  more  thoroughly  for  a  safer  or  more  effective  design, 
and  government  (regulators  and  monitors)  to  oetter  understand  the 
problems  associated  with  this  and  other  unusual  aircraft 
configurations  and  flight  conditions  in  order  to  establish 
specifications,  criteria  and  regulations. 

These  goals  are  seen  to  be  as  follows: 

1.  Upgrade  computational  capability  which  is  intended  to 
support  real-time  flight  simulation  with  graphics,  a 
cockpit,  and  associated  instrumentation  and  controls 
with  force  feedback.  It  is  also  intended  to  support  the 
computat ional  solutions  of  aerodynamic  loads.  This 
includes  the  setup  overall  computing  system:  interfaces; 
D/A,  A/D  converters,  servos  for  stick  force,  graphics 
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drives,  communication,  etc.  and  the  development  of  more 
experience  with  the  system  by  attending  seminars,  and 
using  software  support  services  to  develop  more 
efficient  computational  schemes. 

2.  Further  development  of  the  Aerodynamic  loads  prediction 
techniques  to  include: 

a.  Very  large  angles  of  attack. 

b.  Large  angular  rates 

c.  Arbitrary  wing  planforms. 

d.  Control  surface  effectiveness. 

e.  tdselage  effects  (wing  wake  effect  on  fuselage). 

f.  Propeller  effects. 

g.  Selection,  acquisition  and  implementation  of  a 
suitable  alternative  to  the  present  paneling  code 
for  high  angular  rates. 

h.  Develop  methodology  for  aerodynamic  computations  on¬ 
line  with  the  motion  simulation. 

i.  Extension  to  transport-type  aircraft  planforms. 

*  Note  that  d  &  e  are  required  to  study  spin  recovery. 

3.  Development  of  "Spin  Warning"  Devices  for  alerting 
pilots  of  conditions  which  may  lead  to  spin  entry,  and 
the  development  of  optimal  spin  recovery  techniques 
(i.e.,  with  minimum  loss  of  altitude.) 

5_i.  CONCLUSIONS. 

1.  Tools  for  evaluation  of  aircraft  configurations  and 
their  proneness  to  spin  entry  have  been  developed  which 
show  significant  progress  toward  more  comprehensive 
evaluation  methods. 

2.  For  the  visual  simulation,  more  efficient  computer  codes 
and/or  higher  performance  devices  (now  available)  are 
desirable  for  better  results. 

3.  For  the  aerodynamic  prediction,  code  improvement  and 
more  validation  with  wind  tunnel  or  other  experimental 
results  are  required.  In  addition,  rotational  rates  to 
simulate  aircraft  rotation  and  spin  condition  should  be 
incorporated . 

6_;_  RECOMMENDATIONS . 

The  recommendations  for  further  work  are  outlined  in  the  proposal 
summarized  in  Appendix  C. 
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APPENDIX  A 


On  Estimating  Aircraft  Nonlinear  Rotary 
Derivatives  From  Static  Wind  Tunnel  Data. 

James  R.  Ritte;* 

The  Wichita  State  University 

ABSTRACT 

One  problem  in  mathematical  modeling  of  aircraft  post-stall 
behavior  lies  in  the  difficulty  in  describing  aerodynamic  moments 

due  to  airplane  rates  of  rotation.  In  particular,  the  rolling 

C  C  C  f 

and  yawing  rate  derivatives  7  ,  7p,  Up  and  rij,  are  crucial  to 

accurate  motion  modeling.  These  derivatives  are  not  constant 
values  when  the  linear  assumptions  of  unstalled  flight  are  aban¬ 
doned.  One  method  for  calculating  these  derivatives  was  applied 
to  a  low-wing,  single-engine  general  aviation  airplane  for  which 
sufficient,  high  angle- of- at tack  wind  tunnel  data  is  available. 
This  method  involves  approximations  that  limit  the  validity  of 
the  model  and  produce  unsatisfactory  results  when  used  in  an 
established  spin  modeling  program  and  compared  with  extensive, 
full-scale  flight  test  results.  By  using  a  new  method  and  an 
extended  computational  technique,  a  better  match  with  flight  test 
data  was  obtained.  The  response  to  controls  during  stall  depar¬ 
ture  is  characteristic  of  the  subject  airplane  when  the  new 
derivatives  are  used.  Unlike  the  old  derivatives,  the  new  ones 
are  sensitive  to  variations  of  the  state  space  parameters,  and  in 
a  predictable  manner.  Sensitivity  to  aerodynamio  spanwise 
loading  is  also  evaluated.  Comparisons  are  made  between  the  old 
and  new  derivatives,  and  between  outputs  of  otherwise  identical 
executions  of  the  spin  program. 
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A  List  o£  Symbols 


wing  span  (feet) 
blockage  factor 


center  of  gravity 
rolling  moment  coefficient, 


Q«to  S  b 


rolling  moment  coefficient,  non-rolling 
airplane 

3-d  lift  force  coefficient 


rolling  moment  due  to  rolling  motion 
coefficient,  a  rate  derivative  ("damping-ln- 
roll"). 


C 


? 


P  = 


i(p  b  /  2  V) 


rolling  moment  due  to  yawing  motion 
coefficient,  a  rate  derivative. 


ac, 


'*r  = 


Mr  b  /  2  V) 


3-d  normal  force  coefficient 


yawing  moment  coefficient, 


q«*S  b 


section  normal  force  coefficient  at  station  i 


yawing  moment  due  to  rolling  motion 
coefficient,  a  rate  derivative, 


3(p  b  /  2  V) 


yawing  moment  due  to  yawing  motion 
coefficient,  a  rate  derivative  ("damping-in¬ 
yaw"), 


i(r  b  /  2  V) 


a  distribution  factor  at  station  i 


Symbols  (continued) 

a  general  force  on  the  left  wing;  direction 
is  implied  from  the  context 

a  general  force  on  the  right  wing;  direction 
is  implied  from  the  context 

a  wing  station  designator  variable 

rolling  moment  (foot-pounds) 

lift  (pounds),  left 

length  from  airplane  c.g.  to  the  vertical 
tail  center  of  pressure  (feet) 

slope  of  a  curve;  the  incremental  change  in 
the  ordinate  with  respect  to  the  abscissa  at 
a  point  on  the  curve 

non-dimensional 

roll  rate  (radians  per  second) 

local  dynamic  pressure  ( pounds/ sqr.  foot) 

freestream  dynamic  pressure  (lb. /sqr.  foot) 

yaw  rate  (radians  per  second) 

right 

revolutions  per  minute 

p 

effective  strip  area  at  station  i  (ft*) 

p 

Total  wing  planform  area  (ft*) 

effective  planform  area,  the  sum  of  all  the 
strip  areas,  one  side  only  (ft2) 

time  (seconds) 

trailing  edge  (of  a  wing  or  control  surface) 

body  motion  velocity  with  respect  to  the 
atmosphere,  direction  is  parallel  to  the 
roll  (x-body)  axis  (ft/sec) 

body  motion  velocity  with  respect  to  the 
'.tmospbere,  direction  is  parallel  to  the 
pitch  (y-body)  axis  (ft/sec) 

velocity  (ft/ second) 
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Symbols  (continued) 


Lroll 


*i 


The  product  of  rotation  rate  and  station  span 
arm  to  the  rotation  axis  (ft/second) 

velocity  at  or  very  near  the  airplane 
surface  (ft/ second) 

velocity  in  the  freestreaa  (ft/ second) 

body  motion  velocity  with  respect  to  the 
atmosphere,  direction  is  parallel  to  the  yaw 
axis  (ft/sec) 

distance  from  the  airplane  center  of  mass  to 
a  wing  station;  may  be  in  feet  or  non- 
dimensional 
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reek  Symbols 


®< 

°<A 

°<W 

$ 

A 


an  aerodynamic  angle- of- at  tack,  the  angle 
between  the  longitudinal  axis  and  the  flight 
path  tangent  in  the  vertical  plane  (degrees) 

angle- of- at tack  at  station  1  (degrees) 

angle-of-attack  of  the  airplane 

angle- of-attack  of  the  wing;  differs  from 
angle-of-attack  of  the  airplane  by  the  wing 
incidence  with  respect  to  the  airplane 
longitudinal  axis  (degrees) 

an  aerodynamic  angle-of-attack,  the  angle 
between  the  longitudinal  axis  and  the  flight 
path  heading;  also  called  the  sideslip  angle 
(degrees) 

change  in 


A-6 


Some  Nomenclature 


alpha 

angle-of-attack 

attractor* 

beta 

flat  spin 

limit  cycle* 

pitch 

rate  derivative 

rectangular 

distribution 

rotary  derivative 
Schrenk 

spin 

spiral 

stability 

derivative 

stall 


angle- of- at  tack  in  the  vertical  plane  (seeo<) 

the  angle  between  a  body  axis  and  the  wind 

an  equilibrium  state  that  attracts 
trajectories;  trajectories  need  not  enter  an 
attractor  to  be  affected  by  it;  trajectories 
need  not  interact  in  a  clearly  defined 
fashion  (chaotic  attractor) 

side  slip  angle  (see ^3) 

a  spin  in  which  the  airplane  attitude  is 
relatively  flat  with  respect  to  the  ground; 
alpha  is  generally  greater  than  1  radian 

a  closed  trajectory;  a  helix  in  the  case  of 
the  path  of  a  spinning  or  spiraling  airplane; 
other  aircraft  motions  not  so  easily  defined 

aircraft  attitude  in  the  vertical  plane 

a  derivative  of  a  moment  with  respect  to  a 
rotary  motion 


a  uniform  distribution  of  forces  across  the 
wing  span 

see  rate  derivative 

a  spanwise  distribution  approximation  that 
averages  the  planform  width  distribution  with 
an  elliptical  one 

a  high  angle-of-attack  helical  flight  path 
of  negligible  helix  radius 

any  helical  flight  path  where  the  helix 
radius  is  large,  generally  greater  than  one 
wing  span 


a  partial  derivative  of  a  force  or  moment 
with  respect  to  a  motion 


a  sudden  separation  of  flow  on  the  wing  that 
typically  results  in  a  flight  condition  that 
is  unstable  in  pitch,  roll  and  yaw 


A-7 


Nomenclature  (eont.) 


state  space  a  geometric  model  for  the  set  of  all 

idealized  states  of  an  organism;  the  state 
space  is  filled  with  trajectories  and  is 
called  the  phase  portrait  of  a  dynamical 
system 

steep  spin  a  nose-low  spin  in  which  alpha  is  moderate, 

generally  less  than  1  radian 


trajectory  the  time  rate-of-change  of  a  state  variable 


« 


These  terms  are  part  of  the  language  of  the  geometry  of 
behavior,  discussed  by  Abraham  and  Shaw1'1 
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On  Estimating  Aircraft  Nonlinear  Rotary 
Derivatives  From  Static  Wind  Tunnel  Data 


James  R ,  Ritter 
The  Wichita  State  University 


Introduction 

Most  flight  regimes  permit  linear  modeling  of  the  aero¬ 
dynamic  forces  and  moments  as  functions  of  the  motion  variables. 
This  allows  the  definition  and  use  of  what  are  called  stability 
derivatives,  which  are  the  partial  derivatives  of  forces  or 
moments,  either  in  dimensional  or  non-dimensional  form,  with  re¬ 
spect  to  motion  variables  such  as  velocity,  angle  of  attack, 
rotational  rates  and  so  forth.  Such  derivatives  are  taken  to  be 
constants  for  a  given  flight  regime  in  the  linear  model. 

In  the  approach  to  stall,  and  particularly  in  the  post-stall 
maneuver,  significant  angular  rotation  rates  develop.  These 
rates  cause  linear  approximations  to  be  invalid,  primarily  be¬ 
cause  of  the  aerodynamic  forces  on  the  wing  resulting  from  roll 
and  yaw  rates.  In  stall/spin  maneuvers,  the  partial  derivative 
constants  of  stability  must  be  replaced  by  forces  and  moments 
which  are  functions  of  more  than  one  argument  of  the  motion 
variables. 

IM  Problem 

Aerodynamic  forces  and  moments  can  be  defined  quantitatively 
in  the  wind  tunnel,  but  once  the  airplane  begins  to  rotate  in 
post-stall  departure,  static-flow  measurements  are  not  directly 
representative.  There  are  numerous  schemes  for  evaluating  ro- 
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tating  forces  and  moments,  either  by  rotating  or  shearing  the 
airflow,  or  by  rotating  the  model  in  a  static  flow  field,  and  all 
add  a  great  deal  of  complexity  to  the  management  and  collection 
of  data.  There  has  been  some  hope  that  rotational  aerodynamics 
can  be  derived  from  static  wind  tunnel  data  Methods  used  here¬ 
tofore  are  attempts  to  retain  some  of  the  simplicity  afforded  by 
stability  derivatives,  but  make  them  functions  of  a  single  argu¬ 
ment  rather  than  constanta  Since  the  wing  is  the  major  contri¬ 
butor  to  motion  during  the  stall/spin,  efforts  will  focus  on  an 
improved  method  for  the  four  rate  derivatives  governed  by  the 

c  c  c  c 

wing  aerodynamics,  \p,  tip,  tr  and  nr 

Consider  the  local  angle  of  attack  across  the  span  of  a  wing 
when  significant  roll  and  yaw  rates  are  present  in  the  airplane 
motion  (Figures  1,2,  and  3) ' 
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Figure  2:  Local  Angle  of  Attack  Alpha 
Dominant  Rolling  Shown 


This  means: 

a. )  local  normal  and  axial  forces  vary  across  the  span, 

b. )  a  spanwise  distribution  must  be  assumed,  typically  as 

per  Schrenk  or  rectangular,  and 

c. )  local  dynamic  pressure  varies  across  the  span  due  to 

velocity  increments  due  to  roll  and  yaw. 

Thus  roll  and  yaw  moments  become  functions  of  four  arguments 
which  must  be  determined  by  integrating  across  the  span  of  the 
wing,  i.e.,  =  Cj(«  ,  V,  p,  r  )  and  Cn  =  Cn(«<  ,  V,  p,  r) 

(see  Appendix  B). 
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The  Wvkes  Equations 

One  formulation  used  for  converting  static  aerodynamics  to 
rotating  aerodynamics  has  already  been  made^,  and  the  ensuing 
derivations  are  hereafter  referred  to  as  the  Wykes  equations. 

The  Wykes  derivations  are  the  result  of  a  quasi-static, 
wing-panel  strip  integration  scheme  that  determines  rate  deriva¬ 
tives  as  single-variable  functions  of  airplane  angle  of 
attack  alpha.  Thi3  method  of  derivation  is  accomplished  by  using 
a  simple,  closed-form  approach  to  the  solution,  and  as  a  result 
the  Wykes  equations  are  independent  of  aircraft  motion  to  a  large 
degree. 

The  researcher  who  first  used  the  Wykes  equations  to  inves¬ 
tigate  the  rate  derivatives  of  a  low-wing,  single-engine  general 
aviation  trainer  presented  his  results  in  a  report^)  composed 
chiefly  of  plots  (hereafter  called  Bibrle  derivatives),  very 
basic  descriptions  of  methods,  and  a  list  of  the  Wykes  equations. 
A  clue  to  more  specific  procedures  oomes  from  an  earlier 
report^  in  which  the  Wykes  equations  were  used  to  determine  the 
rate  derivatives  of  a  Navy  trainer.  In  that  report,  sample 
derivations  are  inoluded,  reflecting  for  example  the  use  of  an 
elliptical  lift  distribution  and  a  strip  integration  resolution 
of  five  strips  per  side. 

Tfce  Ritter  Functions 

As  a  simple  and  expedient  means  of  distinguishing  between 
methods  to  obtain  derivatives,  the  rotational  aerodynamic  moments 
developed  in  this  thesis  will  hereafter  be  called  Ritter 
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functions.  The  Ritter  functions  vary  in  accordance  with  the 
premise  that  the  difference  between  wing  panel  moments  is  depend¬ 
ent  on  the  rates  themsalves  and  on  the  velocity  of  the  airplane 
in  a  specific  way.  The  Ritter  functions  depend  on  the  chosen 
state  variables  of  the  system  in  a  non-differentiable  way;  conse¬ 
quently,  these  functions  are  not  themselves  derivatives.  They 
define  moment  coefficients,  replacing  those  that  are  normally 
defined  in  the  spin  equations  by  multiplying  the  rates  with  the 
rate  derivatives.  The  Ritter  technique  requires  only  static  wind 
tunnel  data,  not  rotary  balance  data  or  data  from  other  methods, 
and  these  static  data  are  available  for  the  subject  airplane^. 

Discussion  Concerning  the  Wvkes  and  Ritter  Methods 

Stall  departure  and  spins  involve  significant  angular  rota¬ 
tion  rates,  and  the  primary  influence  of  rotational  motion, 

particularly  in  roll  and  yaw,  comes  from  the  variation  in  local 

r 

angle  of  attack  across  the  span  of  the  wing.  The  variable 
which  in  coefficient  form  represents  rolling  moment  (1)  caused  by 
rolling  motion  (p)  about  the  longitudinal  axis,  is  chosen  to 
compare  the  methods  of  Vykes  and  Ritter. 

r • 

In  the  development  of  aerodynamic  forces  across  the 
span  of,  and  normal  to,  the  planform  of  the  wing  are  considered 
to  act  perpendicular  to  their  respective  span  distances  from  the 
rolling  axis  to  produce  a  net  rolling  moment.  The  continuous 
nature  of  the  force  distribution  is  approximated  by  partitioning 
the  span  of  the  wing  into  strips  and  either  assuming  a  spanwise 
force  distribution  such  as  elliptical  or  rectangular  or 
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Schrenkian,  or  by  using  measured  pressure  data^  from  which 
force  data  can  be  derived  (Appendix  A).  When  the  moments  are 
integrated  by  summation  over  left  and  right  wing  semi-spans, 
a  net  differential  moment  is  the  result.  The  calculation  is 
repeated  over  and  over  for  a  full  range  of  airplane  alphas 
(angles-of-attack)  which  the  aircraft  is  expected  to  en¬ 
counter  during  stall  departure  and  spin. 

In  the  derivation  of  the  equations  of  motion^,  the  rate 
derivatives  were  assumed  to  be  single  variable  functions  of  angle 
of  attack  alpha,  as  per  Wykes.  The  moment  coefficients  and  CQ 
were  reconstituted  from  these  assumed  derivatives  and  their 
respective  angular  rates  during  motion  Integration: 


C2  *  C*n.r.  +  ♦  C*r(«)«r  and 

Cn  a  cnn>r<  «■  Cnp(«)xp  +  Cnr(«x)xr 


For  the  Ritter  method,  the  moments  are  calculated  more  directly. 
In  essence,  local  angle  of  attack  at  each  ith  wing  station  is 
found  from  V,  p,  and  r;  local  dynamic  pressure  is  calculated  for 
each  station;  local  section  normal  and  axial  force  coefficients 
are  found  at  eaofa  station;  force  and  subsequently  moments  are 
found  at  the  ltb  station  and  then  summed  over  both  wing  panels. 
The  equations  are: 


C2  S  cz(i,  V,  p,  r)  .- ZKCHL  y  h 

and  Cn  *  Cn(*,  7,  p,  r)  =.  I^CA  (£)\<  S 7. 


where  R  and  L  refer  to  right  and  left  wing  panels. 
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Vt 

v> 


r 

=  a*  ^°<i)  is  the  nornal  force  coefficient  at 

wing  station  i, 

Vj. 

-yi—  (p,  r)  is  the  local-to-freestream  velocity  ratio, 

®  representative  of  the  dynamic  pressure  ratio, 


y^^  is  the  moment  arm  from  station  i  to  the  roll  axis  of  the 
airplane, 


d i  is  a  distribution  factor  at  station  i  (Schrenk,  elliptical, 
rectangular,  experimental,  etc.), 

e^^  is  the  effective  load  area  for  wing  panel  strip 
representation,  and 

swe  is  the  total  planform  wing  area  Sw  minus  the  area  blocked  by 
the  fuselage,  and  minus  other  areas  not  represented  in  the 
strip  integration  scheme. 

Details  of  the  derivations  may  be  found  in  Appendix  B. 


This  kind  of  multi-variable  representation  means  that 
highly  repetitive  and  often  redundant  calculations  are  required 
for  determining  stability  moments  —  the  elegant  but  less-than- 
accurate  method  of  tfykes  must  be  replaced  with  an  extended 
computer  technique. 


A  Sgaale  SgBEaciaaB 


The  difference  between  Wykes  and  Ritter  can  be  described  in 
the  context  of  the  following  figure  (Figure  4). 


The  above  figure  is  a  plot  of  the  wing  normal  force  coef¬ 
ficient  Cjj  as  it  varies  with  airplane  angle- of- at  tack  o<  .  Point 
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A  occurs  at  example  In  the  Wykes  method,  the  slope  m  at  A 

is  used  to  project  the  Cjj  condition  at  each  wing  partition 
location.  This  method  is  satisfactory  at  points  B  or  C  where  the 
slope  is  more  or  less  constant,  or  anywhere  the  roll  rate  and 
thus  the  o(  variation  is  very  small.  But  in  the  region  around 
point  A,  the  slope  is  not  constant.  The  airflow  becomes 
increasingly  unsteady  on  both  wing  panels,  left  and  right; 
intuitively,  the  resultant  normal  forces  on  each  panel  may  be 
quite  different  at  any  given  instant,  resulting  in  a  large  moment 
one  way  or  the  other,  and  subsequently  a  high  roll  acceleration. 

A  roll  rate  of  150  degrees  per  seoond  (Figure  8)  for  an  airplane 
traveling  at  a  speed  of  120  feet  per  second  (Figure  6)  means  wing 
tip  speeds  of  (tt/  180 )*(  150)x<  1 2.25)  =  32.1  feet  per  second 

(Vron  =  <o  *  v  )  and  an  increased  (or  decreased)  angle  of  attack 
of  approximately  15  degrees: 

A<k  ?  arotan  (  Vroll/V ^  )  a  arctan  (  32.1  /  120  )  =  15  degr. 
The  Wykes  derivative  thus  fails  to  be  consistently  representative 
sinoe  the  difference  between  wing  tip  alphas  can  be  as  great  as 
30  degrees.  Dynamic  pressure  also  varies  in  a  way  that  compounds 
the  problem  for  positive  alpha,  positive  roll  rate  (Appendix  B). 

Wykes'  method  is  concluded  by  taking  the  derivative  of  roll 
moment  with  respect  to  roll  rate  and  plotting  the  results.  In 
contrast,  the  Ritter  method  is  ideally  concluded  at  the  deriva¬ 
tion  of  the  net  moment.  There  is  no  way  to  plot  Ritter  functions 
that  is  consistent  with  the  premise  that  each  moment  coefficient 
be  evaluated  at  each  instant  in  time  with  respect  to  all  relevant 
state  variables,  but  the  Ritter  functions  can  be  evaluated  and 
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compared  with  the  Bihrle  derivatives  if  a  set  of  fixed  conditions 
is  chosen.  A  baseline  of  variables  is  selected,  in  which 
velocity  is  120  feet  per  second  and  roll  rate  is  one  radian  per 
second.  A  pseudo-derivative  curve  (C^/p,  CQ/r,  etc.)  is  then 
produced  that  is  representative  of  the  selected  states.  The 
curves  of  Ritter  and  Bihrle  "derivatives"  are  plotted  and  com¬ 
pared  in  Figure  5. 

The  same  method  used  to  calculate  the  [  function  is  used 
to  calculate  c£r,  Cnp,  and  cnr  functions.  The  four  derivatives 
of  the  Wykes  method  are  shown  in  Appendix  C  and  compared  with  the 
specific  parametric  evaluations  of  the  Ritter  functions  as  in 
Figure  5.  (The  Wykes  and  Ritter  methods  might  be  expected  to 
compare  favorably  when  the  roll  rate  is  small;  a  roll  rate  of  .1 
rad/ second  is  evaluated  and  the  comparison  is  as  expected.) 

One  unknown  condition  that  required  investigation  was  spanwise 
distribution.  Integrated  pressure  distributions^  and  two  as¬ 
sumed  distributions,  rectangular  and  Schrank,  are  compared  in 
Appendix  D.  There  are  some  small  differences  in  two  of  the 
derivatives  when  pressure  data  is  compared  with  the  assumed 
methods.  These  differences  are  attributed  to  inconsistencies 
between  force  and  pressure  data,  and  to  the  varying  flow  in  the 
alpha  range  investigated  (11.2  to  41.5  degrees).  There  is  no 
noticeable  difference  between  reotangular  and  Schrenk  distribu¬ 
tions  when  considering  a  baseline  set  of  state  variables 
(V  =  120  ft/sec  and  p  and  r  =  1  rad/second). 
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The  Evaluation 

The  complex  behavior  of  an  aircraft  in  stall  or  in  spin  can 
be  described  symbolically  as  the  solution  of  a  set  of  first-order 
differential  equations  in  which  both  the  kinematic  and  aerody¬ 
namic  terms  involve  nonlinearities.  In  solution,  no  limit  cycle 
can  be  readily  described;  nevertheless,  actual  spin  behavior  is 
mode-like,  and  the  attractor  or  attractors  (many  airplanes  have 
more  than  one)  can  be  most  easily  evaluated  through  numerical 
integration  and  plotting. 

A  great  many  full  scale  flight  test  time  histories  exist  for 
the  subject  aircraft  from  NASA  spin  research^.  This  means  an 
actual  result  is  available  to  prove  or  disprove  a  model.  A 
six-degree-of-freedom  model  of  the  equations  of  motion^ 10^  has 
been  used  for  five  years  at  Nichita  State  University  in  spin 
research^1  These  equations  are  without  restriction  in  the 
limits  of  motion,  and  contain  fully  nonlinear  kinematic  terms. 

The  Ritter  functions  provide  a  complementary  set  of  compatible, 
nonlinear  aerodynamic  terms. 

A  program  was  written  to  evaluate  the  Ritter  method.  It  was 
used  independently  to  generate  the  specific,  single- variable 
derivatives  in  Figure  5  and  in  Appendices  C  and  D.  An  evaluation 
of  spin  behavior  that  is  consistent  with  the  premises  concerning 
the  state- dependent  nature  of  the  derivatives  as  functions  is 
made  by  incorporating  this  program  into  the  spin  program  as  a 
subroutine,  thus  replacing  piecewise  linear  fits  of  the  Bihrle 
derivatives.  The  spin  solutions  using  the  Bihrle  derivatives  and 
the  Ritter  functions  are  then  compared  with  each  other  and  with 
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flight  teat  data  (Figures  6,  7,  8,  and  9). 

A  flight  test  spin  was  chosen  which  reflects  a  varied,  full- 
range  application  of  elevator  and  rudder.  These  control  profiles 
have  been  modeled  in  the  spin  program  and  plots  of  three  key 
trajectories,  o( ,  p,  and  r,  are  compared  on  the  following  pages. 
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gure  6:  Flight  Test  of  a  Low-Wing,  Single-Engine 
General  Aviation  Airplane  in  Spin: 

Speed,  Power,  Elevator  and  Rudder  Profiles 
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Figure  7:  Flight  Test  vs.  Bihrle  and  Ritter  Spins: 
Airplane  Angle-of-Attack  Alpha 
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The  Bihrle  derivatives,  based  on  the  Wykes  equations,  were 
used  in  a  flight  simulator  program  as  a  means  of  qualifying  spin 
behavior.  The  response  to  controls  had  been  shown  to  be  very 
sensitive  in  early  studies;  nevertheless,  the  variable  deriva¬ 
tives  were  effective  in  producing  a  spin.  The  Ritter  functions 
have  produced  a  more  realistic  departure  and  response  to  controls 
up  until  the  point  where  power  was  added  during  flight  test 
recovery.  Engine  thrust  has  not  been  modeled  in  the  spin 
program;  hence,  a  Limit  of  Validity  has  been  established  in 
Figures  6  through  9  defining  that  point  in  time. 

The  improved  behavior  using  Ritter  functions  is  ascribed 
at  least  in  part  to  the  expanded  concept  used  in  defining  the 
rate  derivatives,  that  the  moments  CQ  and  Cg  are  functions  of 
,  p,  r,  and  V  and  can  be  defined  numerically  and  in  a  consistent 
manner  from  within  the  research  spin  program. 
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Appendix  A:  Pressure  Data  and  Wind  Tunnel  Force  Data 

Researchers  at  NASA  collected  wing  pressure  data  for  the 
full-scale  test  airplane  airplane  angles-of-attack  11.2,  21.6, 
31.9,  and  41.5  degrees^ A  pressure  integration  program  was 
written  to  evaluate  the  normal  and  axial  forces  at  the  six 
apanwise  locations  selected  for  pressure  testing. 

The  pressure  tests  are  substantive  along  the  chord  and  are 
used  to  evaluate  absolute  pressure  forces  in  the  region  of  each 
spanwise  station.  A  unit  width  of  one  inch  is  considered  at  a 
point  of  measurement,  surface  length  segment  boundaries  are  shared 
between  pressure  ports,  and  the  product  of  width  segment 
and  length  segment  is  the  local  area  at  each  port. 

A  B-spline  curve  fit  is  used  on  the  upper  and  lower  surfaces 
to  evaluate  the  slope  at  all  points  respective  to  the  chord 
reference  line.  The  pressure  forces  are  first  evaluated  perpen¬ 
dicular  to  the  surface.  The  components  of  each  force  normal  to 
the  longitudinal  axis  of  the  airplane  are  collected  along  the 
section  line  and  summed  to  yield  a  normal  pressure  force  on  the 
wing  at  that  wing  station.  Similarly,  axial  components  are 
collected  and  summed.  The  result  is  a  set  of  six  normal  and  six 
axial  forces  across  each  wing  panel,  left  and  right  sides.  The 
profile  of  the  forces  in  each  direction  is  a  representation  of 
the  force  distribution. 

In  order  that  these  distributions  be  in  some  way  compatible 
with  distributions,  Schrenk  or  otherwise,  that  must  be 
assumed  outside  the  pressure  tested  alpha  range,  the  pressure 
distributions  are  normalized  to  a  value  of  one.  The  general 
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shape  of  the  pressure  force  distribution  is  a  guide  to  what 
assumed  distribution  might  be  appropriate. 

Sohrenlc's  approximation  was  initially  used  for  the 
normal  distribution  above  and  below  the  pressure  range. 

Below  11.2  degrees,  is  approximately  equal  to  C^.  Above 
41.5  degrees,  Schrenk  is  used  because  the  pressure  force 
distribution  at  41.5  degrees  is  roughly  similar. 

A  rectangular  distribution  is  used  for  axial  distribution. 
Below  11.2  degrees,  2-d  friction  drag  is  more  significant  than 
pressure  drag,  and  axial  forces  are  similar  by  the  small  angle 
relationship.  At  41.5  degrees  the  distribution  is  flat  with  only 
a  slight  increase  (negative  forward)  at  the  wing  tip. 

Figure  A1  shows  the  normal  and  axial  pressure  force 
distributions  after  normalization.  The  lift  and  drag 
distributions  are  included  for  reference  in  Figure  A2.  Actual 
forces  are  not  shown  since  they  are  not  compatible  with  the  wind 
tunnel  force  coefficient  data  used  throughout  the  angle- of -at tack 
range  of  investigation  (-8  to  90  degrees) ^ . 

The  wind  tunnel  force  data  is  shown  in  Figure  A3.  The 
normal  force  represents  a  combination  of  a  full-scale  wind  tunnel 
evaluation  from  -8  to  41  degrees  angle-of -attack  and  a  1/7  scale 
model  evaluation  beyond  41  degrees.  Fuselage-only  data  is 
subtracted  from  wing-fuselage  data  and  the  resulting  wing  data  is 
corrected  for  fuselage  blockage  using  a  blockage  factor 
b.f .  =  1/(1  -  (blocked  area  /  unblocked  area)  ).  The 

axial  data  is  treated  the  same  way. 
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Figure  A2:  Lift  and  Drag  Pressure  Force  Distributions 
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Appendix  B:  The  Ritter  Functions 


1.)  The  Derivation  of  Cj(o<,  7,  p)  (traditionally 
"Damping  in  Roll") 


ci s  ct(a<  ’v>  p)  = 

R  and  L  refer  to  right  and  left  wing  semi-spans. 
r 

is  the  normal  force  coefficient  at  station  i, 

"^T  is  the  local-to-freestream  velocity  ratio, 

y^  is  the  moment  arm  from  station  i  to  the  roll  axis  of  the 

airplane, 

dj^  is  the  distribution  factor  at  station  i, 

e ^  is  the  effective  load  area  for  strip  representation,  and 

we  is  the  total  planform  wing  area  minus  the  area  blocked 
by  the  fuselage  and  other  areas  not  represented  in  the  strip 
integration  scheme. 

These  variables  are  discussed  below  in  the  order  presented. 


is  the  local  normal  force  at  station  i.  Since  CN  is 
a  function  of  airplane  angle- of -attack,  can  be  looked  up  if 

is  known.  Cjj  of  the  wing  is  determined  from  three- 
dimensional  wind  tunnel  data  (Appendix  A)  and  is  originally 
distributed  across  the  span  in  a  uniform  way.  The  actual 
spanwise  strength  distribution  is  adjusted  by  another  factor,  d^ 
Station  angle-of-attack  is  determined  from  the  airplane 
angle -of- at tack  plus  the  change  in  angle-of-attack  due  to  the 
rolling  velocity,  V1  roll  :  =  ®<A  +  A*c 
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Figure  B1 


The  total  local  velocity  is  affected  at  each  point  on  the  wing. 


(Wing  Incidence  is  neglected  since  force  data  was  evaluated  at 
airplane  alpha. )  vdWar*/  \\*6 
For  the  left  wing  ^ 


^4  =  0X^0 


the  angle-of- 
attack  decreases. 


The  non- rolling  alpha  is  defined  by  the  state  variable 
alpha.  The  rolling  A *1  is  defined  by  the  state  variables  Vm  and 
p  as  follows: 

Right  wing:  I 


A*/ 
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Using  the  law  of  sines, 


-Acxt  ) 


f  V 


v« 


and  through  trig,  identities,  A«</  ='fck« 


-j 


Vo,  +  P^t  <>"  <*  J 


Similarly  for  the  left  wing: 


=  t<+vC 


n 


;  C6S  o( 


L  *'*<*.  j 


(^sp^O 


Figure  B5 


It  can  be  shown  that  these  equations  are  valid  for  positive  and 
negative  combinations  of  roll  and  yaw.  It  only  necessary  to 
remember  that  for  the  left  wing,  o( ^  s 

The  most  effective  way  to  define  Cjj  vs.  a(  within  the 
program  is  from  a  table  of  values.  This  means  digitizing  of 
the  appropriate  curves  from  reference  5.  This  was  done  by 
selecting  a  few  critical  points  and  fitting  a  cubic  spline  to 
them.  Graphic  comparisons  are  made  to  insure  a  smooth  and 
satisfactory  fit  (Appendix  A).  The  ensuing  table  of  cublcs  is 
initiated  into  an  array  in  a  program  used  to  evaluate  the 
derivatives.  As  conditions  are  encountered,  the  table  can  be 
read  and  interpolated  as  required. 

Since  it  is  the  local  dynamic  pressure  that  accounts  for  the 
actual  forces  (  Fn  =  CQ  q  S  )  the  ratio  of  local  pressure  to 
freestream  pressure  is  required:  ^ 

't  f— 
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so  it  is  the  ratio  of 


For  the  conditions  investigated  ^ L 
velocities  that  is  all  that  needs  to  be  considered.  For  veloc¬ 
ity  ratio  expressions,  it  is  better  to  use  the  law  of  cosines, 
thus  avoiding  singularities  possible  with  the  law  of  sines.  The 
following  expression  has  been  shown  to  be  v— 11J  for  all  comoina- 
tions  of  roll  rate  and  angle-of-attaclc. 

/  Vk.  \  f  pu«  fp>4.  +  2  \Jot  4lVt°<  ) 

(-£ )  -  1  +  r  v - 

For  the  left  wing,  the  second  sign  in  the  expression  is  changed. 
The  velocity  ratio  for  the  left  wing,  positive  alpha,  positive 
roll,  is  seen  to  be  less  than  one. 

yi  (non-dimensional) 

The  pressure  tests  of  reference  6  were  oonduoted  at  semi- 
span  locations  2y/b  a  .38,  .44,  .63,  .78,  .85,  and  .92.  These 
values  are  adopted  for  this  analysis.  Actual  lengths  are  used 
in  the  scalar  approach  of  the  preceding  derivations. 

d^  (non-dimensional) 

These  factors  represent  the  spanwise  distribution  across  the 
effective  use  area  of  the  wing  (next  paragraph).  All  distribu¬ 
tions  have  been  normalized  to  have  a  mean  value  of  one 
(Appendix  A.) 

Self  Effective  Load  Area  for  Integration 

Forces  on  the  wing  oause  moments  about  the  axis  of  rotation. 
To  get  a  fair  calculation,  at  least  the  outer  one-half  of  each 
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wing  panel  must  be  considered.  Rather  than  interpolating  a 
uniformly  distributed  array  of  span  locations  to  be  evaluated  at 
uniform  areas,  tha  selected  pressure  regions  are  used  with  areas 
defined  within  shared  boundaries. 


Figure  B6 


CD  ®  ®  ®  © 


Table  B1 

Column  0  values  are  normalized  wing  stations  ~ 

Column  0  values  are  half  the  distance  between  stations. 

Column  0)  values  are  shared  boundary  stations. 

Column  (7)  values  are  actual  areas  associated  with  station  i 

locations,  ^e^.  ^we  is  the  sum  of  column  0}  values 
in  Table  B1.  Swe  =  29.65  ft2. 

Column  (T)  values  are  normalized  area  ratios  ^e^/^w@. 
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The  Rolling  Moment 


The  total  moment  is  evaluated  as  follows: 


Figure  B7 

When  p  is  positive,  FR  will  produce  a  negative  moment  and  a 
positive  moment.  When  the  airplane  is  flying  much  below  stall, 

FR  will  be  greater  than  FL  and  the  net  moment  due  to  rolling  will 
be  negative  (stabilizing). 

The  product  of  force  and  moment  arm  coefficients,  with  the 
necessary  adjustment  coefficients  for  dynamic  pressure,  area  and 
distribution,  give  the  total  moment  coefficient  evaluated  for 
each  wing  panel,  left  and  right.  The  3um  of  left  and  right 
total  moments  gives  a  net  rolling  moment  due  to  rolling 
coefficient: 


V*’  V,  P)  =  - C„, (-£)•  J. ^ ^ (£);■* 

r 

An  artificial  derivative  * i p  can  be  constructed  for  the 

purpose  of  comparing  Ritter  functions  with  Bihrle  derivatives,  i.e., 

c  c 

Jp  =  ?/p  where  p  is  non-zero.  When  used  as  a  subroutine 

in  the  spin  integration  program,  it  is  only  the  moment  itself 
that  is  required,  so  conditions  of  zero  roll  rate  do  not  present 
the  problem  of  indefinable  equations.  When  the  subroutine  is 
used  independently  as  a  research  tool  the  roll  conditions  are 
controlled  to  avoid  the  problem. 

The  following  page  describes  the  different  configurations 
for  angle-of-attack  and  rolling  motion. 
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Right  wing,  positive  alpha,  positive  roll  rate: 


Figure  B8 


V. 


alpha  becomes 
more  positive 


Right  wing,  negative  alpha,  positive  roll  rate: 


Figure  B9 


alpha  becomes 
less  negative 


Right  wing,  positive  alpha,  negative  roll 


alpha  becomes 
less  positive 


Figure  BIO 

Right  wing,  negative  alpha,  negative  roll  rate; 


•V 


Figure  B1 1 


alpha  becomes 
more  negative 
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2.)  The  Derivation  of  C j(-<,  V,  r)  (traditionally  cir,  Roll 
Moment  Due  to  Yaw  Motion) 

The  correct  equations  are  derived  from  the  following  dia¬ 
grams  in  a  manner  similar  to  the  derivations  in  Part  1.)  of  this 
Appendix. 

Right  wing,  positive  alpha,  positive  yaw  rate: 


more  positive 

Figure  B12 

Right  wing,  negative  alpha,  positive  yaw  rate: 


Right  wing,  positive  alpha,  negative  yaw  rate: 


/  less  positive 


Figure  B14 
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Right  wing,  negative  alpha,  negative  yaw  rate: 


W^ateffcieu^  N 


Figure  B15 


alpha  becomes 
less  negative 


■[a  Vi 

r  u  5'»» °<  l 

L  —  lA  Coy 

A<s<L  - 

faw" 

U  5iv\o(  1 

L  V<*  ul  Cj>Sa{  J 

c*L 

=  ®< 

♦  A«<<. 

=  <?< 

-  A  *Kc  (Left  wing) 

Only  the  first  expression  above  is  required  to  evaluate  all 
combinations  of  r  and  e<  for  the  right  wing. 

Cj  S  C7(  «(,  V,  r)  *  -T  f . •  J +  S'  /  .  •  f 


T 
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3.)  The  Derivation  of  C  ( <* ,  v,  p)  (traditionally  ^n_t  Yaw 
Moment  Due  to  Roll  Motion) 

cnp  differs  from  ^7 p  and  ''ip  in  that  axial  force  data  is 
used  to  calculate  the  net  moment.  CA  of.  the  wing  is  determined 
by  subtracting  CA  body  wind  tunnel  data  from  CA  data. 

The  resulting  curve  is  seen  to  be  negative  (thrusting)  between  5 
and  21  degrees  alpha  (Appendix  A). 

Cn  =  Cn^  * '  V»  p)  s 
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4.)  The  Derivation  of  Cn(  e( ,  V,  r)  (traditionally  unr, 
"Damping  in  Yaw") 


r  c 

For  both  np  and  nr,  caution  must  be  used  in  interpreting 


data.  Aa  aeen  in  Appendix  A,  drag  ia  alwaya  poaitive  (rearward) 
but  ttte  axial  force  diatribution  ia  quite  variable  in  polarity. 

cn  =  cn(af'  v’  r)  =  5^ 


A-42 


Appendix  C:  Bihrle  Derivatives  and  Ritter  Functions  Compared 


The  following  charts  reflect  direct  evaluation  of  the  Wykes 
equations  and  conditional  evaluation  of  the  Ritter  equations. 
Baseline  conditions  for  evaluation  of  the  Ritter  equations  are, 
velocity  is  120  feet  per  second  and  roll  and  yaw  rates  are  both 
one  radian  per  second. 

The  Bihrle  derivatives  for  the  subject  airplane  are  repro¬ 
duced  from  the  report^  ^  and  shown  in  Figure  Cl.  The  Ritter 
baseline  evaluations  are  displayed  in  Figure  C2. 

Because  the  Ritter  equations  are  dependent  on  more  than  one 
flight  parameter,  they  can  only  be  conditionally  displayed.  The 
following  pages  show  the  variations  in  the  functions  resulting 
from  different  velocities  (Figure  C3)  and  from  different  roll  and 
yaw  rates  (Figure  C4). 
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The  Bihrle  Derivatives 


Non -dimens  Cot  f. 


This  chart  is  constructed  on  half-aegree  aloha  increments 


Ritter  derivatives,  V/120/.  Beta/0/,  P/1/,  R/l/ 


Figure  C2a:  Baseline  Ritter  Functions 
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-  C  1  r 

-  Cnr 

This  chart  is  constructed  on  half-degree  alpha  increments. 


Ritter  derivatives,  V/i20/,  Beta/0/,  P/1/,  R/l/ 


Figure  C2b:  Baseline  Ritter  Functions,  Reduced  Scale 
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Figure  C3a:  The  Ritter  Functions,  V  =  80  ft/sec 
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This  chart  is  Constructed  on  half-degree  alpha  increments. 
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Figure  C3b:  The  Ritter  Functions,  V  =  1 60  ft/sec 
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-10-5  0  5  10  15  20  25  30  35  40  45  50  55  b0  b5  70  75  80  85  90 

Airplane  Angle  of  Attack  Alpha  (degrees) 

Ritter  derivatives.  V/120/.  Beta/0/,  P/.  1/,  R/.  1/ 

Figure  C4a:  The  Ritter  Functions,  p  and  r  *  .1  rad/sec 


Cnp  -  C  I  r 


This  chart  is  constructed  on  half-degree  alpha  increments. 
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Figure  C4b:  The  Ritter  Functions,  p  and  r  =  .5  rad/sec 


Non  d  i 
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Figure  C4c:  The  Ritter  Functions,  p  and  r  =  1.8  rad/sec 
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This  chart  15  constructed  on  half-degree  alpha  increments. 


Airplane  Angle  of  Attack  Alpha  (degrees) 
R.tter  derivatives,  V/120/,  Beta/0/,  P/2.5/.  R/2. 5/ 


Figure  C4d:  The  Ritter  Functions,  p  and  r  =  2.5  rad/sec 


Appendix  D: 


Sensitivity  to  Snanwise  Force  Distribution 

A  conditional  set  of  Ritter  "derivatives"  (V  =  120  ft/sec, 
p  and  r  =  one  rad/sec)  were  evaluated  for  different  spanwise 
distributions. 

The  original  research  distribution  (Schrenk  for  normal 
forces,  rectangular  for  axial)  is  used  to  derive  the  rate  deriva¬ 
tives  shown  in  Figure  C2  of  Appendix  C..  The  pressure  distribu¬ 
tion  described  in  Appendix  A  is  applied  in  computing  the  deriva¬ 
tives  shown  in  Figure  D1.  The  Schrenk  distribution  is  applied 
for  all  derivatives  in  Figure  D2  and  the  rectangular  distribution 
is  applied  for  all  derivatives  in  Figure  D3. 
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Non 


Pressure  data  spanwise  distribution 
all  derivatives 
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Figure  Dla:  Baseline  Ritter  Functions 
Pressure  Force  Distribution 
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This  chart  is  constructed  on  half-degree  alpha  increments. 
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Figure  Dlb:  Baseline  Ritter  Functions,  Reduced  Scale 
Pressure  Force  Distribution 


Non-dimens  Coef. 
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Rectangular  distribution,  a 
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Figure  D3a:  Baseline  Ritter  Functions 
Rectangular  Distribution 
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This  chart  is  constructed  on  half-degree  alpha  increments. 


Reciangu  i-ar  distribution,  all  derivatives 
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Figure  D3b:  Baseline  Ritter  Functions,  Reduced  Scale 
Rectangular  Distribution 
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EFFECTS  AT  HIGH  ANGLES  OF  ATTACK 
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Abstract 

A  method  is  presented  for  the  estimation  of  the 
effects  of  wing  planform  modifications  on  its  aero¬ 
dynamic  forces  using  the  popular  panel  methods. 
The  objective  of  such  estimation  is  to  assist  in 
the  preliminary  concept  design  of  spin  resistant 
aircraft.  The  pressure  distributions  are  obtained 
by  applying  changes  to  the  boundary  conditions  at 
the  panel  control  points  to  simulate  separated 
boundary  layer  flow.  The  changes  in  normal  velo¬ 
city  vectors  are  computed  based  on  2-dimensional 
airfoil  data  beyond  the  onset  of  separation. 
Results  Tor  angles  of  attack  up  to  27°  are  in 
agreement  with  the  experiment.  Another  run  at  a 
32°  angle  of  attack  shows  a  deviat.on  from  experi¬ 
mental  data,  however,  the  error  is  the  same  for  the 
baseline  and  the  planform  modification  considered. 

nomenclature 


b  wing  span 

C]_(y)  local  lift  coefficient  at  a  wing  section 
at  a  span  location  y 

C,  lift  curve  slope  in  the  linear  range 

C*  pressure  coefficient 

L,  »;,N  aerodynamic  rolling,  pitching  and  yawing 

moments 

n  normal  vector  to  surface 

P,Q,R  roll,  pitch  and  yaw  rates 

V  velocity  vector 

x  chordwise  distance  from  leading  edge 

y  distance  from  the  plane  of  symmetry 

y'  non-dimensional  y  =  2y/b 

angle  of  attack  of  the  entire  wing 
a(y)  local  angle  of  attack  at  a  span  location  y 

p  bound  filament  vortex  strength 

a ,  source  strength  for  panel  1 

Subscripts 


1  or  J  panel  or  section  number 
n  normal  component 

1.1.  limit  of  linear  range  of  C ^ 

sep  separation  point  on  airfoil 

o  induced  due  to  singularities 

»  free  stream 


Introduction 

The  stall/spin  characteristics  of  general  avia¬ 
tion  aircraft  have  been  the  subject  of  renewed 
interest  in  recent  years.  Flight  safety  consider¬ 
ations  have  led  to  a  thorough  flight  test  and  other 
experimental  programs  at  NASA-Langiey  and  a  number 
of  universities. 

Spin  entry  results  at  hi#i  angles  of  attack  due 
to  the  resulting  adverse  rates  of  change  of  aero¬ 
dynamic  moments  with  respect  to  angular  velocities 
of  the  aircraft.  These  moments  are  caused  mainly 
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by  surface  pressure  distributions  on  the  wing,  but 
other  surfaces  trailing  it  (fuselage  tailcone  and 
stabilizer  surfaces)  are  also  involved.  Even 
though  their  contribution  i3  very  Important,  these 
other  surface  effects  are  caused  by  the  wing  wake 
which  becomes  unsymmetric  upon  spin  entry.  There¬ 
fore,  much  of  the  spin  research  concentrates  on 
wing  character  at  high  angles  of  attack. 

In  designing  spln-resi3tant  aircraft  configur¬ 
ations,  a  preliminary  step  is  to  examine  the  lift 
coefficient  vs.  angle  of  attack  curve  beyond  stall. 
A  sharp  drop  In  lift  with  increasing  a  indicates 
the  potential  for  unfavorable  rolling  moment  rate 
of  change  with  roll  rate  (dL/dp).  A  strip  (or 
preferably  surface)  integration  of  the  moments  of 
normal  and  axial  forces  along  the  wing  span  is 
necessary  for  a  more  conclusive  evaluation  of  wing 
character,  3ince  the  effective  angle  of  attack 
varies  with  spanwlse  location'.^. 

For  a  conf iguration  to  be  spin  resistant,  or  at 
least  to  enter  the  spin  gently,  the  wing  planform 
must  be  such  that  the  post-stall  behavior  of  it3 
outboard  panels  is  mild  (the  outboard  panels  are 
more  critical  due  to  their  larger  moment  arm). 
This  is  illustrated  in  Figure  1. 


Figure  1.  Effect  of  local  m  on  lift  and  hence 
rolling  moment. 


In  order  to  achieve  this  design  objective,  an 
outboard  panel  leading  edge  cuff  (extension)  was 
developed  and  tested  at  HAEA-Lar.gley 3.  It  was  also 
adopted  for  at  least  one  recent  general  aviation 
design.  The  abrupt  change  in  chord  length  and 
leading  edge  camber  causes  a  delay  in  the  stall  of 
the  outboard  panel  which  provides  the  desired  spin 
character.  Reference  3  3hows  oil  flow  visualiza- 
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tlonn  of  an  aircraft  wing  with  and  without  the 
leading  edge  device  near  and  beyond  C^max.  The  oil 
flow  photographs  clearly  show  the  outboard  panel 
remaining  attached  for  the  modified  wing. 

The  objective  of  the  present  work  is  to  enable 
computational  prediction  of  this  behavior  in  order 
to  provide  preliminary  designers  with  a  capability 
to  evaluate  planform  modifications. 

The  method,  which  is  in  its  early  stages  of 
dev elo pmen t ,  is  based  on  a  special  application  of 
potential  flow  (panel)  methods  to  simulate  sep¬ 
arated  now.  The  idea  is  to  artificially  generate 
a  boundary  layer  effect  by  altering  the  surface 
boundary  condition  in  the  solution  of  Laplace's 
equation.  Panel  methods  solve  this  equation  by 
enforcing  boundary  conditions  of  zero  flow  across 
solid  boundaries.  In  order  to  simulate  the  bound¬ 
ary  layer,  a  nonzero  normal  velocity  distribution 
on  the  separated  portions  of  the  wing  surface  is 
specified  as  boundary  condition. 

The  method  presented  is  simple,  requires  minor 
modifications  to  existing  potential  now  codes  with 
which  many  aerodynamic  configuration  designers  are 
familiar.  It  provides  good  estimates  of  the  change 
in  character  due  to  planform  modifications. 


Following  the  success  of  the  leading  edge  cuff 
design  and  the  results  of  its  oil  flow  aerodynamic 
evaluation,  a  detailed  study  [4]  was  conducted  at 
the  Wichita  State  University  wind  tunnel  facility 
to  understand  the  mechanism  of  delaying  separation 
on  the  outboard  panel  of  such  a  wing. 

In  this  study,  velocity  and  pressure  data  were 
obtained  with  a  five  tube  pressure  sensing  probe. 
Measurements  in  the  flow  field  above  and  behind  the 
wing  upper  surface  were  taken  so  that:  a)  vortl- 
city  could  be  detected,  and  b)  the  separated  wake 
could  be  determined. 


Figure  2.  Cuffed  wing  configuration. 


This  investigation  led  to  the  following  conclusions: 

].  The  sudden  change  in  chord  at  the  notch  leads 
to  a  change  in  bound  vorticity  (see  figure  2). 

2.  In  addition,  especially  at  high  a,  the  peak 
negative  pressure  near  the  leading  edge  is  much 
lower  outboard  of  the  notch  due  to  the  droop. 

3.  The  combined  effect  of  the  above  leads  to  the 
shedding  of  a  strong  vortex  from  the  leading  edge. 


4.  This  vortex  interacts  with  the  tip  vortex  and 
is  pulled  toward  the  tip,  introducing  energy  into 
the  boundary  layer  and  reducing  its  thickness 
considerably. 


Figure  3.  Vortex  Interaction. 


5.  The  flow  does  not  really  remain  attached  on 
the  outboard  panel,  but  is  leas  gravely  separated 
with  the  introduction  of  the  cuff. 


Comoutatlona 


In  the  present  work,  the  separated  flow  is 
modeled  only  by  modifying  the  boundary  condition  of 
the  potential  flow  panels.  No  attempt  was  made  to 
account  for  the  interaction  of  the  two  trailing 
vortices.  This  would  require  an  iterative  (relax¬ 
ation)  scheme.  Also,  the  vortex  sheet  is  left  on 
the  surface  of  the  wing.  A  more  realistic  model 
would  allow  for  its  displacement  away  form  the 
surface  in  the  separated  regions.  The  normal 
velocities  used  to  raodiry  the  panel  boundary  con¬ 
ditions  to  simulate  a  boundary  layer  are  computed 
for  chordwise  2-D  airfoil  sections.  In  thav  model, 
a  single  vortex  filament  at  the  quarter  chord  is 
used  rather  than  a  distributed  vortex  sheet  and 
crossflows  are  not  accounted  for.  These  simpli¬ 
fying  assumptions  lead  to  some  estimation  error, 
but  adequately  allow  prediction  of  the  Incremental 
effects  of  small  changes  in  planform. 


Description 


the  Method 


The  computation  is  performed  in  the  following 
■u^uence: 

a)  A  3-D  potential  flow  analysis  is  performed  on 
the  wing  at  the  desired  angle  of  attack  with  no 
simulated  separation.  The  objective  ‘s  to  obtain 
an  initial  CL(y )  distribution  by  strip  integration 
of  the  local  pressure  coefficients  (Figure  4a). 

b)  The  flow  is  assumed  strictly  chordwise.  For 
each  section  i  (located  at  y^),  a  local  angle  of 
attack  a(yj)  Is  obtained  which  corresponds  to 
C1  Cy i )  by  solving  the  equation  of  the  linear  por¬ 
tion  of  the  lift  curve  (the  dashed  line  in  Figure 
4b),  where  the  slope  is  . 

This  local  a(y)  1®  *n  initial  estimate  of  the 
finite  span  effect  of  the  wing  at  the  local  section. 

c)  If  the  flow  Is  separated  at  the  section  in 
question  (i.e.  if  the  local  a  exceeds  the  linear 
limit  o1-1_  then  c^yi)  Is  read  from  Figure  4b 
(solid  line)  as  well  as  x3ep  and  Cp3ep  (figure  4c), 
These  three  values  must  be  available  from  2-D  alr- 
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foil  data  (experimental  or  computational). 


(y ' ) 


Figure- 4a.  Potential  flow  spanvise  lift 
distribution 


Cf. 


Cj^2  *  Required  for  a(yj.) 


Figure  4b.  Airfoil  lift  curve 


Figure  4c.  Required  separation  data  for  local 
airfoil  section 


d)  The  assumption  is  made  that  the  pressure 
coefficient  remains  constant  downstream  of  the 
separation  point.  The  problem  now  is  to  "design" 
the  circulation  and  the  boundary  conditions  for  the 
airfoil  section,  namely  Vn^  for  each  control  panel 
on  the  separated  flow  surface  so  that: 

1.  The  lift  coefficient  is  equal  to  the  finite 
wing  local  c^(y^)  and, 


2.  The  pressure  coefficient  at  all  panel  con¬ 
trol  points  downstream  of  the  separation  point  is 
equal  to  Cpsep  corresponding  to  a  j,. 

This  is  performed  iteratively  by  using  5  function 
minimization  technique  (steepest  decent).  The 
problem  is  then: 

minimize  f  =  K[C1-C1(y1) ]2  +  j  [CpJ-Cpsep]2 

*j>xsep 

where  and  Cpj  are  functions  of  the  design 
variables  T  ,  V n j .  They  are  computed  for  each 
variation  of  these  design  variables. 

The  function  f  represents  the  difference  between 
the  desired  separated  pressure  distributions  and 
lift  coefficient  and  the  ones  resulting  from  the 
computations.  and  Cpj  are  computed  using  a  2-D 
panel  program  which  features  a  single  vortex  of 
strength  T  at  the  aerodynamic  center  and  constant 
source  panels.  The  panel  source  strengths  and  the 
circulation  r  are  computed  in  the  manner  of  panel 
methods  by  imposing  a  3et  of  boundary  conditions  of 
specified  normal  velocity  Vn  at  each  panel  control 
point,  plus  the  Kutta  condition  at  the  trailing 
edge  (Figure  4d).  The  solution  is  reached  when  f 
becomes  small  or  after  a  finite  number  of  iter¬ 
ations. 

Tangential  Flow 


Figure  4d.  2-D  Airfoil  panel  method  analysis 

e)  The  normal  velocitie*  computed  for  each  span 
section  are  then  used  to  modify  the  boundary  condi¬ 
tions  in  the  3-D  potential  flow  program  used.  The 
aerodynamic  influence  coefficient  matrix  [AIC]  i3 
unchanged.  The  vector  o  is  obtained  by  solving 
the  equation 

[AIC]  [o  ]  »  [b] 

where  a  is  the  unknown  singularity  vector  (source 
and  doublet  strengths)  and  b  is  the  vector  con¬ 
taining  the  desired  normal  velocity  components  at 
all  control  points.  The  normal  velocity  vector  Vni 
at  panel  1  control  point  is 

Vni  =  (V.  «■  V0  )  •  n 

Vni  =  0  for  attached  flow  and  i  0  for  separated 
flow.  The  specific  values  used  are  those  obtained 
from  steps  d. 

f)  The  modified  boundary  condition  run  results  in 
a  re-computed  spanwise  lift  distribution  which 
accounts  for  separation.  The  process  may  be 
repeated  starting  from  step  b,  except  that  the 
solid  line  in  Figure  4b  could  now  be  U3ed  instead 
of  the  dashed  line,  since  "better"  values  of  sec¬ 
tion  lift  are  available. 
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Results  and  Discussion 
First  Phase  of  Development 

The  Initial  development  of  the  present  work  was 
based  on  experimental  results  published  in  refer¬ 
ence  5.  There,  a  full  scale  wind  tunnel  test  of  a 
general  aviation  aircraft  (complete  configuration) 
was  conducted.  Several  leading  edge  oonf igurations 
were  added  and  tested  also.  These  results  com¬ 
prised  pressure  distributions  at  various  span 
stations  and  lift  coefficient  curves  up  to  a=ao°. 
These  were  used  to  test  the  validity  of  the  present 
method  in  very  general  terms. 

A  panel  model  for  the  wing  alone  was  used  in  the 
coraput  .tions,  the  fuselage  was  omitted  for  the  sake 
of  expediency.  Some  discrepancy  was  expected  with 
the  wind  tunnel  results  for  that  reason.  At  that 
point  in  the  development,  only  one  iteration  was 
performed,  i.e.,  the  pure  potential  flow  analysis 
was  done,  followed  by  the  computation  of  the  normal 
velocities  required  to  simulate  the  boundary  layer 
thickness,  based  on  the  potential  flow  lift  distri¬ 
bution  Ci(y),  and  finally  a  recomputation  of  the 
wing  lift  was  done. 

The  computation  was  performed  for  the  basic  wing 
and  one  leading  edge  modification  which  has  been 
tested.  The  curves  in  figure  5  are  reproduced  from 
reference  5.  Also  shown  In  figure  5  are  the 
results  of  the  computations  shown  with  solid 
circles  (basic  wing)  and  squares  (modified  leading 
edge). 

In  these  computations,  2-D  airfoil  data  were 
required.  Since  there  were  no  detailed  2-D  data 
available  for  the  tested  airfoil  section,  the  pres¬ 
sure  data  plots  published  in  reference  5  were  used. 
For  the  basic  airfoil,  the  chordwlse  Cp  distribu¬ 
tion  of  the  section  closest  to  the  plane  of  sym¬ 
metry  was  used.  For  the  cuffed  airfoil,  the  sec¬ 
tion  at  the  midpoint  of  the  wing  outboard  panel 
with  the  leading  edge  euff  was  used.  These  were 
considered  close  approximations  of  airfoil  data 
even  though  they  were  taken  from  a  finite  wing. 
From  the  Cp  plots,  available  for  a  range  of  angles 
of  attack,  the  location  of  the  separation  point  and 
the  pressure  coefficient  from  there  to  the  trailing 
edge  were  plotted  as  functions  of  a,  and  used  for 
the  estimation  of  the  "separated"  boundary  condi¬ 
tions. 

A  look  at  the  total  configuration  lift  coeffi¬ 
cient  for  three  angles  of  attack  of  21.6°,  27°  and 
31.9°  (figure  5)  Indicates  that  up  to  27°,  che 
computational  model  agrees  with  the  experiment, 
whereas  at  31.9°  it  did  not,  but  the  deviation  is 
nearly  the  same  for  the  two  planforos  which  were 
compared. 

The  results  at  ct  =  3 q .9°  are  valid  In  the  sense 
of  evaluating  the  incremental  effect  of  small  plan- 
form  modifications,  a  capability  which  is  useful  in 
designing  such  changes,  with  the  intent  of  allevi¬ 
ating  the  sharp  drop  in  lift  after  stall. 

In  figure  6,  the  computed  pattern  of  the  start 
of  separation  for  various  angles  of  attack  is 
shown.  On  the  outer  panel,  the  separation  occurs 
further  downstream  when  the  leading  edge  cuff  is 
included.  These  plots  are  consistent  with  the  oil 
flow  photographs  of  reference  3. 

These  preliminary  results,  despite  all  the 
short-comings  In  the  data  used,  were  encouraging 
enough  to  pursue  this  work  further. 


Figure  5.  Comparison  of  computed  and  experimental 
lift  coefficients  of  test  wing. 


Figure  6.  Computed  stall  pattern  of  basic  and 
cuffed  wings. 


Second  Phase  of  Development 

New  wind  tunnel  tests  were  started  at  Wichita 
State  University  using  a  simple  rectangular  win; 
with  a  NACA  24023  airfoil  section  and  an  aspec- 
ratio  of  13.3.  In  this  test,  pressure  data  wer‘ 
collected  using  a  Strip-a-Tube  belt. 

These  tests  are  still  underway  at  the  time  o:' 
this  writing,  therefore  only  the  analysis  of  tha 
basic  wing  is  performed.  The  objective  is  to  tee¬ 
the  convergence  of  the  iterative  scheme.  The  larg- 
aspect  ratio  and  the  fact  that  there  is  no  fuselag- 
helps  determine  the  separation  point  and  the  2-1 
lift  coefficient  more  reliably.  Plots  of  X.ep  an: 
C„,.„  for  various  angles  of  attack  are  depicted  ir 
figures  7  and  8. 

The  ehordwise  pressure  was  integrated  to  produc- 
the  experimental  local  lift  coefficient.  The  C 
vs.  a  curves  for  various  y  locations  are  consider¬ 
ably  different  In  magnitude,  but  the  shape  i: 
similar.  Airfoil  data  exists  for  this  airfoil  ir 
reference  6  and  was  used  because  of  the  measureraer. 
error  existing  1  the  data  due  to  the  relatlv 
thickness  of  the  .  .rlp-a-Tube  belt. 
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and  lta  corresponding  a.  See  figure  9.  When  the 
Ci  values  are  within  the  feasible  range,  then  the 
corresponding  a's  are  taken  from  the  curve. 


Figure  7.  Location  of  the  separation  point. 


Figure  9.  Intermediate  lift  curves  for  the 
computation. 


e)  By  performing  4  iterations,  the  values  of 
appear  to  be  converging  toward  the  experimental 
data.  The  f.nal  results  of  the  computation  are 
compared  with  experiment  in  figure  10.  The  history 
of  the  convergence  is  shown  in  figure  11. 


Figure  3.  Upper  surface  pressure  coefficient  in 
the  separated  region. 

The  computational  method  was  tested  as  follows: 

a)  The  panel  method  model  was  constructed  for  the 
wing  using  ten  rows  of  chordwise  panels. 

b)  A  pure  potential  flow  (Vn  s  0)  run  was  per¬ 
formed  for  a  wing  angle  of  attack  of  30°.  The  lift 
at  the  section  nearest  the  plane  of  symmetry  was 
considered  the  2-0  lift  coefficient  (value  z  3.07) 
at  a  =  30°.  At  a=0,  CL  =  0.10.  A  straight  line 
interpolant  was  used  to  obtain  the  local  angle  of 
attack  a(y)  for  determining  x3ep  and  Cp3ep. 

c)  The  7n  vector  was  computed  for  each  spanwise 
section  and  another  potential  flow  run  was  made 
with  the  new  boundary  conditions.  These  V„  are 
based  on  the  potential  flow  lift  distribution  and 
hence  account  for  downwa3h. 

d)  Some  of  the  resulting  section  lift  coefficient 
values  cause  a  problem  in  finding  the  corresponding 
a  since  is  limited  in  range,  except  for  the 
linear  portion  of  the  curve.  For  example,  a  of 
0.6  or  less  would  imply  a  low  local  a  (attached 
flow)  even  though  it  is  known  that  the  flow  is 
highly  separated  at  this  point.  On  the  other  hand, 
a  very  large  C|  is  not  realizable  and  no  correspon¬ 
ding  u  exists.  Iteration  of  the  procedure  thus 
necessitates  a  different  approach  for  determining 
the  3(y)  for  the  next  iteration. 

The  values  of  C^  fall  outside  the  feasible  range 
frequently  because  the  ct's  are  over  predicted  or 
under  predicted  by  the  potential  flow  run,  depend¬ 
ing  on  the  extent  of  the  separation.  Therefore, 
using  the  average  Cj_  of  the  previous  and  current 
estimates  is  appropriate  to  help  accelerate  conver¬ 
gence.  In  addition,  the  largest  value  of  =<long 
the  span  is  used  along  with  the  wing  angle  of 
attack  to  serve  as  one  point  on  a  2-piece  linear 
Interpolant  for  o.  The  other  two  points  required 
are  the  azO  and  its  corresponding  Clt  and  Cimax 
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Figure  10.  Computed  and  experimental  lift 
distribution. 


Conclusions 

The  method  presented  here  is  being  developed  as 
a  part  of  a  program  dealing  with  stall/spin  safety. 
It  is  Intended  to  provide  the  preliminary  designer 
with  a  quick  way  of  determining  effects  of  wing 
planform  variations.  Adding  angular  velocities  to 
more  clearly  evaluate  the  aerodynamic  moments  of 
the  wing  can  be  done  by  further  modifying  the 
boundary  conditions. 
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5.  Newsom,  W.A.  Jr,  et  al,  "Efrects  of  Wing 
Leading  Edge  Modifications  on  a  Full-Scale,  Low- 
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Wing  Sections,"  Dover  Publications,  New  York,  1959. 


Figure  11.  Convergence  of  the  computed  solution. 


The  method's  known  weaknesses  are  those  due  to 
some  simplifications  adopted  for  expediency.  For 
example,  in  the  estimation  of  the  normal  velocity 
components,  a  single  lifting  line  approach  was 
used,  rather  than  a  distributed  vorticlty.  This 
led  to  3ome  oscillations  of  pressure  distributions 
near  the  trailing  edge,  and  a  large  pressure  peak 
at  the  leading  edge.  Further,  the  displacement  of 
the  vortex  sheet  away  from  the  upper  surface  and  a 
better  application  of  the  Kutta  condition  should 
improve  the  results.  More  runs  to  compare  with 
existing  data  are  required  to  refine  and  validate 
the  method. 

Its  strengths  are  its  siapli city,  and  the  fact 
that  it  uses  the  reliable  and  familiar  potential 
flow  codes.  The  computing  requirements  are  reason¬ 
able.  So  far,  the  comparisons  with  experimental 
data  Involved  many  approximations.  It  is  desirable 
to  conduct  more  testing  with  pressure  taps  for  a 
variety  of  planforms.  The  results  presented  show 
much  potential  despite  these  weakness. 
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APrENDIX  C 


- WICHITA  STATE  UNIVERSITY  HZ 

ZZZ  INSTITUTE  FOR  AVIATION  RESEARCH 


Project  Title:  Stail/Spin  Flight  Simulation,  Phase  II:  Spin  Resistant  Design  and  Spin  Entry  Detection 

Principal  Investigator( a):  M.G.  Nagati 

Problem: 

•  Lack  of  a  design  methodology  and  tools  for  the  design  of  spin  resistant  aircraft  (non-linear,  high 
angle  of  attack  and  angular  rate  aerodynamic  analyses). 

•  Need  for  a  “spin  warning-’  device.  Development  requires  the  knowledge  of  parameters’  behavior  at 
onset  of  spin. 

•  Realistic  evaluation  of  warning  systems  and  spin-resistant  designs  requires  flight  simulator  with 
special  characteristics. 

Objective: 

•  Develop  aerodynamic  analysis  methods  for  wing,  fuselage,  empennage  in  the  presence  of  high  angular 
rates,  separation  and  propeller  wash. 

•  Identify  parameters  most  dominant  in  spin  entry  and  study  behavior  under  a  variety  of  conditions 
Accomplishments  to  date: 

•  Good  simulation  of  aerodynamic  forces  and  aircraft  dynamic  response  has  been  established. 

•  Good  prediction  or  stall  pattern  for  isolated  wing  with  and  without  leading  edge  extensions  has  been 
accomplished. 

•  Survey  has  been  done  of  propeller  wake  prediction  methods  for  twin  engine  aircraft. 

•  Graphics  for  flight  simulation  have  been  developed. 

Brief  work  statement: 

1.  Complete  simulator:  install  Learjet  cockpit  and  add  visual  effect  (computer  graphics),  checkout  by 
test  pilots. 

2.  Perform  runs  to  identify  spin  entry  parameters;  develop  requirements  for  “spin  warning-’  device. 

3.  Select,  acquire,  and  implement  computational  aerodynamic  codes  for  arbitrary  shapes  with  separated 
flows.  Modify  codes  to  allow  for  angular  rates. 

4.  Continue  development  of  modified  panel  method  to  obtain  effects  of  leading  edge  shape  on  wing 
stall,  and  extension  of  technique  to  more  arbitrary  platforms. 

5.  Study  and  incorporate  effects  of  fuselage  and  wing  wakes  on  aerodynamic  forces  during  spin. 

6.  Incorporate  propeller  induced  flow  effects. 

T.  Publish  results. 

Estimated  costs:  $140,000  per  year. 
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